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Abstract 

We consider the boundary crossing problem for time-homogeneous diffusions 
and general curvilinear boundaries. Bounds are derived for the approximation error 
of the one-sided (upper) boundary crossing probability when replacing the original 
boundary by a different one. In doing so we establish the existence of the first- 
passage time density and provide an upper bound for this function. In the case of 
processes with diffusion interval equal to M this is extended to a lower bound, as 
well as bounds for the first crossing time of a lower boundary. An extension to some 
time-inhomogeneous diffusions is given. These results are illustrated by numerical 
examples. 
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1 Introduction 

Calculating the probability a given diffusion process will cross a given one- or two-sided 
boundary in a finite time interval is an important problem for which no closed-form 
solution is known except for a few special cases. Such boundary crossing probabilities 
arise in many applications, including finance (pricing of barrier options) and sequential 
statistical analysis. Since no closed form exact solution is immediately forthcoming in the 
general case, finding approximate solutions, together with approximation rates, becomes 
an important alternative approach. A possible pathway here is to replace the original 
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boundary with a close one, for which obtaining boundary crossing probabihties is feasible, 
and then to obtain a bound for the error caused by using the approximating boundary. 

This approach was taken in [5], where an error bound was established for boundary 
crossing probabilities for the Brownian motion. It was shown there that under mild 
regularity conditions, the difference between the probabilities doesn't exceed a multiple 
of the uniform distance between the original and approximating boundaries, with the 
coefficient being an explicit function of the Lipschitz coefficient of the boundary. An earlier 
(asymptotic) bound, again in the Brownian motion case, was obtained in [13], but under 
some additional superfluous conditions. In [16], the authors dealt with a special class of 
diffusion processes which can be expressed as piecewise monotone (not necessarily one- 
to-one) transformations of the standard Brownian motion. For such processes, boundary 
crossing problems could be reduced to similar ones for the Brownian motion process. 
Further references can be found in [13] and [T6] . 

Much work has also been done in the area of calculating the density of the first crossing 
time. Explicit formulae can be obtained for a limited number of specific pairs of diffusions 
and boundaries, such as in [9], [7|, [10] (see also references therein), where the existence of 
the density is usually assumed. Existence of the first crossing density was established in 
[8J under certain conditions on the diffusion and boundary. For ffat boundaries and under 
certain smoothness conditions on the diffusion and drift coefficients (which, for example, 
do not apply to the Bessel process), establishes the existence and smoothness of the 
first-passage density. 

In this paper we extend the work of ^ to general diffusion processes. As an auxiliary 
result (which is also of independent interest), we establish the existence of, and bounds 
for, the density of the first passage time under mild conditions on the boundary and 
diffusion. Most of this work is focused on the time-homogeneous case, however we outline 
how some time-inhomogeneous processes can be treated as well. Section [2] introduces the 
processes involved and some other notation. Section [3] proves the existence of a density 
for the first passage time, and provides upper and lower bounds for the density. This is 
used in Section H] to bound the difference between the crossing probabilities of two close 
boundaries. Section [5] provides some numerical examples. 

2 Setup and Notation 

First we consider time-homogeneous diffusions; a possible generalisation is discussed at 
the end of Section [SI We begin with a diffusion Ut governed by the stochastic differential 
equation (SDE) 

dUt = u{Ut)dt + a{Ut)dWt, (1) 
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where Wt is a Brownian motion and cr{y) is different iable and non-zero inside the diffusion 
interval (that is, the the smallest interval ICR such that Xt E I almost surely). As is 
well-known, one can transform the process to one with unit diffusion coefficient. This is 
achieved by (see e.g. [13], p. 161) defining 



for some t/q from the diffusion interval of Ut, and then considering Xj := F{Ut). By Ito's 
formula, this process will have unit diffusion coefficient and a drift coefficient ^{y) given 
by the composition 



Any boundaries being considered must also be transformed accordingly. From here on we 
work with the transformed diffusion process Xf governed by the SDE 



Conditions mentioned throughout refer to the transformed process X^ and its drift coef- 
ficient ^. 

Without loss of generality, we consider boundary crossing probabilities over a time interval 
[0, 1]. The general case of an interval [0, T] follows by a change of scale. 

An 'upper' boundary for the process X^ is a function g(t) such that g{0) > x, while 
a 'lower' boundary satisfies g{0) < x. When considering a two-sided boundary crossing 
problem with an upper boundary (7+ and lower boundary g^, we assume that g+{t) > g-{t) 
for t G [0,1]. 

We will consider the following two cases only: 

[A] The diffusion interval of Xt is the whole real line R. 

[B] The diffusion interval of Xt is R+ = [0, 00). 

Where the specific interval is not important, the lower end-point (either —00 or 0) will 
be denoted by I. When considering boundary-crossing problems for diffusions on M4., we 
assume the boundaries do not touch zero in the time interval [0, 1]. The results extend to 
diffusions with other diffusion intervals with one finite endpoint by employing appropriate 
transforms. 

We denote by P^; probabilities conditional on the process in question {Xt or some other 
process, which will be obvious from the context) starting at x. Where no subscript is 





dXt = fi{Xt)dt + dWt, Xq = X. 



(2) 
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present, either conditioning is mentioned explicitly or the process is assumed to start 
from zero (the latter will only apply to the Brownian motion process). 

For the diffusion Xt we will need a 'reference' diffusion with certain characteristics. These 
are: being able to bound the probability that this reference diffusion crosses a linear 
boundary when its final value is known; the process has unit diffusion coefficient; and 
the process has the same diffusion interval as X^. If X^ has a diffusion interval of M, 
we use the Brownian motion as the reference process. For diffusions on M+ we use the 
Bessel process of an arbitrary integer-valued dimension d > 3. Note that a further reason 
for using Bessel processes is that they are of independent interest, having applications in 
many areas. 

Recall the definition of the Bessel process Rt of dimension d = 1,2,.... This process gives 
the Euclidean distance from the origin of the d-dimensional Brownian motion, that is, 

where the W^^^ are independent standard Brownian motions, i = 1, . . . ,d. As is well 
known (see e.g. [H], p. 445), Rt satisfies the SDE 

dRt = ^^dt + dWt. (3) 

Note that for non-integer values of d the Bessel process of order d is defined using the 
above SDE. The process has a transition density function pfi{t,y, z) given by 

PR{t. y.z) = z ' r ie-(^^+^^)/^% (^) , (4) 

where t] = d/2 — 1 and Tn{z) is the modified Bessel function of the first kind. For further 
information, see pr4J. 

Finally, by $ we denote the standard normal distribution function. 

3 The Existence and Bounds of the First Passage 
Time Densities 

This section proves the existence of, and gives an upper bound for, the density of the first 
passage time 

T := inf{t >0:Xt> git)}. 

The diffusion Xt, with Xq = x, is subject to some mild conditions and the upper boundary 
g{t) satisfies Lipschitz condition. For diffusions on R the result is then extended to lower 
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boundaries, and further results give lower bounds for the first hitting time density for 
both upper and lower boundaries. 

Let ^{t) denote the maximum of the boundary function up to time t: 

git) := max 

()<s<t 

We assume Xt has a transition density p{t,y,z) and define functions q{t,y,z), G{y) and 
M{t) as follows, according to the diffusion interval of Xf. 

[A] If the diffusion interval of Xt is M, then we denote by q{t, y, z) the transition 
density of the Brownian motion, and define, for some fixed y^ G M, 

ry 

G{y) ■■= I fi{z)dz, 

Jyo (5) 
M{t) := inf (f,'{y) + f^'{y)). 

-x<y<g{t) 

[B] If the diffusion interval of Xf is M+, then we denote by q{t,y,z) the transition 
density of a (i- dimensional Bessel process for some fixed integer d > 3, and define, 
for some fixed yo >0, 



G{y) := r (f,{z) -^]dz, 



vo 



2z 



id-l){d-3) , 



(6) 



M{t) := inf fi'iy) - ^ ^ + f,\y) 

To simplify the statement of our main theorem, we set 

{X in case [Al 
—X m case [BJ. 

Theorem 3.1. Let Xt be a diffusion process satisfying ([2]) with n{y) differentiahle inside 
the diffusion interval. Let g{t) he a function on [0, 1] such that g{f)) > x and, for some 

-K-h< g{t + h) - g{t) < K^h, <t <t + h <1. (8) 

Then the first passage time r of the process Xt of the boundary g{t), t G (0, 1), has a 
density Pr{t) , satisfying 

p^{t) < Big,t), (9) 

where B{g,t) is given by 

B{g,t) := -{g{t) + K-t -x)q{t,x, g{t))e^^3{t))--G{xy{t/2)M(t) ^ ^^q^ 

t 
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Remark 3.2. i) For diffusions on M+, the choice of the integer d is arbitrary, subject 
to d > 3, and can be used to optimise the bound. 

ii) The bound is sharp: for constant fi and a hnear boundary, equahty holds in ((91). 

iii) In terms of the original process f/* in ([T]), the condition of differentiability of fi{y) 
requires ^{z) to be differentiable and a{z) to be twice different iable. 

Proof. For a fixed t and small /i > 0, we can write, conditioning on the value of Xf, 

r9(t) 

f,{t e{t,t + h)) = / f,{t e {t,t + h)\Xt = z)F^{Xt e dz) 

Je 

g(t)~h^/* pgit) 

+ =:h + h. (11) 

Clearly, by the Markov property 

P.(r G (t, t + h)\Xt = z) =P,. ( sup (X, - g{s)) <QXt = z 

\0<s<t 



X P, sup (X, - g{s)) >OXt = z]. (12) 

\t<s<t+h / 

We proceed by bounding the first factor on the right hand side for the region g{t) — h^/'^ < 
z < g{t) (it is not required for the remaining interval, as will be seen later). The process in 
question is now a 'diffusion bridge', or 'pinned diffusion': we are interested in the behavior 
of the original process X^ over < s < t conditional on its value at time t. We will use 
the following lemma. 

Lemma 3.3. Let P^ denote the law of Xs starting at Xq = x and pinned by Xt = z. If 
Ts := o [Xu '■ u < s), then for any A E J-'t 

p[t, X, z) 

where 1b is the indicator of the event B and denotes expectation with respect to the 
probability defined as follows: 

(i) if the diffusion interval of Xg is M, then denotes the law of the Brownian motion 
Wg on [0,t], starting at Wq = x and conditional on Wt = z, while 

N{t) := [\fi'{X^)+fi\X^))du; 
Jo 
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(ii) if the diffusion interval of Xs is [0,oo), then denotes the law of the Bessel 
processes Rg of dimension d > 3 on [0,t], starting at Rq = x and conditional on 
Rt = z, while 

Nit) := 1^ (/.'(XJ - ^^'^^^-^^ + AXu?j du. 



Proof. In case [A], the result was obtained in [2j, so we only need to consider case [B]. 
For the latter, we will use a modification of the argument leading to (14) in 

Let C = C([0, t], M"*") be the space of all continuous positive functions on [0,t], which we 
will interpret as a sample space for the processes in question, and on C define Vs{uj) = Us, 
u & C, J^s = aiVu, u < s). On the space (C, J-'t), in addition to the laws P^^ and P^, define 

Qx '■= the law of Rg, 

Ql := the law of Rg, pinned hj Rt = z, 

where Rs is the (i-dimensional Bessel process starting at Ro = x (see ([3])). Hence under 

Px 

dV, = fi{Vs)ds + dW„ 

while under Q^^ 

dVs = - — -ds + dWs, 
2Vs 

for a Q Brownian motion Wg- Let denote expectation with respect to Q^;. If 

d-w 1 r / .... d- 1^ ^ 



Cs := exp <; / ( /.(K) - - 2 7^ l^-"^^") " ^ ' '^'^ 

exp W ( /^(K) - ^) (rfK - ^-^du 
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then by Girsanov's theorem, for all A G we have 

P,(A) = [CtMv^A]] ■ 

Under Q^^ we have by Ito's formula, using ([6]), that 



so 



Ct = exp <^ G{Vt) - G{Vo) - - / l^'iVu 



d-1 



d-l 

"2Kr 



du 



exp <^ G{Vt) - G(Vo: 



{d- l){d-d,) 



4^2 



+ ^^\Vu)]du\. 



Using Lemma 3.1 of we have 



Inserting ( |T3l) into ( |T^ gives 

exp <! - G(\/o; 



completing the proof of the lemma. 



Next we apply Lemma 13.31 to our boundary crossing problem. 
Lemma 3.4. (i) In case [A], 



P,( sup (X, - < 

0<s<f 



X, = z 



P{t,X,Z) \0<s<t 



w. 



(ii) in case [B], 



F4 sup (X, - < 

0<s<t 



Xt = z 

q{t, X, z) G{z)~G{x)-(t/2)M{t) 



p{t,x,z) 



( sup <0 

0<s<t 



where M{t) was defined in ([5]), ([6]) according to the diffusion interval of Xg. 
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The assertions of the lemma immediately follow from that of Lemma 13.31 with A = 
{supo<s<t — g{s)) < O}. Note that M{t) gives a lower bound on the integrand in 
N(t) for all paths in A. 



In order to apply the above lemma, we introduce the linear boundary 

gtis):= git) + K~{t~s), < s < 1. 
Due to assumption (IHl), we have gt{s) > g{s) for < s < t, and hence 



(15) 



( sup {Ys - g{s)) < 

•.0<s<t 



Yt = z] sup (Ys - gt{s)) < 

0<s<t 



Y 



for any process Yg. We now need to bound the boundary crossing probabilities for Brow- 
nian bridges and pinned Bessel processes for linear boundaries. In the case of Brownian 
bridges, an explicit formula is well-known (see e.g. [Sj, pp. 64-67): 



sup {W, - gtis)) < 



0<s<t 



Wt = z 



l-e 



-2h{x)b2{z) 



(16) 



where bi{x) = {gt{0) — x) and &2(^) = t ^ {gt{0) — K t — z). The pinned Bessel process 
satisfies the following inequality: 

Lemma 3.5. For a d- dimensional Bessel process Rg, d = 1,2,..., and x G {0,g{0)), 
zeiO,git)), 



P,. ( sup {Rs -gt{s))<0 

•.0<s<t 



Rt = z] < 1 - e 



-2bi{-x)b2{z) 



with bi{x) and 62 (-z) defined above. 

Remark 3.6. If the conditioned diffusion (or 'Bessel bridge') Vg, < s < t, in Lemma [3.51 
(whose distribution coincides with the conditional law of Rg, < s < t, given Rt = z) 
followed, by analogy with the case of the Brownian bridge, a SDE of the form 



dr.. 



z -r, d- 1\ , 

+ — — 1 ds + dWs 



t - s 



2r. 



(as stated in pQ), we could obtain a sharper bound than the one from Lemma [3.5[ Un- 
fortunately, it appears that the Bessel bridge process actually satisfies the following SDE 
(which can easily be derived from Theorem 5.8 and (4.dl)-(4.d3) in [12]): 



drg 



^ z 1^ A* ^ — rs \ ds -\ — ds + dWs 



t - s 



2r. 



(17) 
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(cf. (jlj)). Alternatively, this SDE could also be derived using an explicit formula for the 
transition density of the Bessel bridge: 



P:, {Rs e dy\Ru = v,Rt = z) = y- 



t-u e-^e-tr^^X.(^)X.(S) 



s — u){t — s) 



g 2(t — u) 



u < s < t, which is easily obtainable from (1.0.6) of [I], p. 373. With (fTTl) we are unable 
to improve the above result. 

Proof. Let Wg denote the d-dimensional Brownian motion {Ws^\Ws'^\ . . . ,Ws'^'^^ with 
initial value Wq = x = (x, 0, . . . , 0), so we can assume without loss of generality that 
R = IIVFJL s > 0. Then 



sup iRs-gtis))<0 

0<s<t 



Rf = z 



sup {\\Ws\\-gt{s))<0 

KO<S<t 



Pa; ( sup {\\Ws\\-9t{s)) <0 
0<s<t 



Z: \\Z\\ = z 

X P.-r iWtedz 



W, 



< sup Fy( snp i\\Ws\\-9tis))<0Wt = z]. 
z,y:\\z\\=z,\\y\\=x \0<s<t J 

Due to the rotational symmetry of the (i-dimensional Brownian motion, we may rotate 
our coordinates so that 2; = (z, 0, . . . , 0). This gives 



Ri = z 



PJ sup {R,-gt{s))<{) 

\0<s<t 

< sup Py( sup {\\W,\\-gt{s))<^ 
y:\\y\\=x \o<s<t 

< sup ¥y I sup (iVj^) - gti-s)) < 

y:\y\<x ' \0<s<t 

= P_, f sup <0 

\0<s<t 

Applying ( |T6i) concludes the proof. 



W, 



(1) 



.,0) 



Wi'^ = z 



□ 
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With X, M(t) and G defined in ©-(CD and for g{t) - h^/^ <z< g{t) we tfius fiave 



P,( sup (X, -(7(s)) <0 

Q<s<t 



Xt = Z] < ^(^' ^' ^) ^G(^)-G(a:)-a/2)M(f) Q _ g-2fei(x)62(^)) 

- p(t,x,z) 



wliere 



A/,(5(,t) := sup 

g{t)-h^/i<w<g{t) 



< sup q(t,x,w)e 

g{t)-hy*<w<g{t) 

9{t) - z 



G{w)-G{x)-{t/2)M{t) 2&i(a:)&2(^) 

p{t,x,z) 



p{t,x,z)' 

g(t,x,u;)e«('")-^(^)-(*/2)*^W-(^7,(0) 



— a; 



Now return to (fT^ and consider tlie second factor, i.e. the probability that the process 
Xg crosses the boundary g{s) during {t,t + h) given Xt = z. Again we bound this 
probabihty by the corresponding hnear boundary crossing probabihty, with the hnear 
boundary defined in (fTSi) . Due to assumption (IHD, we have gt{s) < g{s) for t < s < 1, and 
hence 



Xt = z]. 



sup (X, - g{s)) >0 Xt = z] sup (X, - gt{s)) > 

\t<s<t+h J \t<s<t+h 

For z G {i,g{t) — h^^^], if the process is to cross the hne gt{s), it must first hit the level 
7 := g{t) — h^/^ some time during + h) (we assume without loss of generality that 
< h^^^^). Thus for these z 



P ( sup (X, - gt{s)) >OXt = z] < sup P sup (X, - gt{s)) > 

yt<s<t+h / t<t'<t+h \t'<s<t+h 



Xf = 7 



<P( sup Xs-gt,h>0 

^t<s<t+h 



Xt = l 



= P^ sup X^ > gt^h 

\0<s<h 

where gt,h ■= gt{t + h) = g{t)-K~h. Now let = rh{X) := inf{s > : X^ > gt^h}- Then 
for an arbitrary a G {i, 7) we have 



P^ sup X, > gt,h ] =^^[rh<K inf X^ > a ] + [ u < h, inf X, < a . 

\0<s<h J V 0<s<r;, J ^ V 0<s<r;, ' 



(19) 



We bound the last probabilities by considering the associated probabilities for processes 
which dominate, or are dominated by (as appropriate) X^ almost surely. For the first 
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term on the right-hand side, consider a process such that Zg > Xs a.s. prior to the 
time t' when Xg leaves the interval {a,g{t)). Namely, we put 

dZs = fids + dWs, s > 0, with jl := sup ii{y) < oo 

due to the continuity of /i inside the diffusion interval. Then clearly 
¥JThiX)<h, inf Xs>a 

V 0<s<rh(X) 



< Th{Z) <h, inf Zs>a] < {Th{Z) < h) 
= P ( sup {Ws + fis) > gt,h - 7 

\0<s<h 

= $ (M^/^ - {g,,, - 7)/.-^/^) + e2^(^*.-^)$ (-^^2 - (^,. - 7)/.-^/^) 
= 1 - $ (/3-(7)) + e^'^^S'^'^-^) [1 - $ 



where the second last equality is a well-known relation for the Brownian motion (see for 
example (1.1.4) on p. 250 in [1]) and f3±{z) := {gt^h — z)h~^^'^ ± fih^/"^. For the second term 
on the right-hand side in (fTOl) . we choose a process Yg such that Yg < Xg a.s. for s < r' . 
Let 

dYg = fids + dWg, s > 0, with fi := inf /u(?/). 

Then clearly 

F^lrhiX) <h, inf X, < a 

V 0<s<rh(X) 



< P^ Th(X) <h, inf Yg < a] { inf < a 

0<s<r;,(X) / ' yO<s<h 



P inf (us + Wg) <-f-a 

\0<s<h 

l-$(A_(7)) + e''^(^-°ni-'^(^+(7))]. 
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where \±{z) := (z — a)h ± fih^/'^. Combining these results gives the bound 



P( sup {Xs-gt{s))>Q 

t<s<t+h 



Xt = z 



< 1 - $ (/3_(7)) + e^^^^'-'^-^^ [1 - $ (/3+(7))] 
+ 1 - $ (A_(7)) + e2^(^-") [1 - $ (A+(7))] 



27r/?-(7) 
1 1 

/27r A_(7)' 



27r/?+(7) 



(20) 



using Mill's inequality. 



Now choose a := g{t) — /i^/^. Then for the first term on the right-hand side of ( |20l) we 
have 

/?_(7) = - K-h - {git) - h'l'))h-'l' - i,h'l' = /.-^^[l - {K- + fc)h'/'], 
and so clearly 

1 e-^^i(7) = 



as h ^ 0. Similarly, as A±(7) = h '^^^{1 + o(l)), other terms are also o{h). Thus for Ji 
in (fTTl) we have 

Ji = o(/i), (21) 
and hence the integral Ji does not contribute to the density of r. 

Next we consider z e [g{t) — /i^/^, (^(t)], using essentially the same argument but without 
conditioning on starting at the point 7, which also eliminates the need to reduce the 
problem to a level crossing. Returning to the original boundary g{t), this gives 

pf sup (X, - g{s)) >OXt = z) <!-<!> {x^^)) + e2('^+^")(^'W-^) [l - $ (^+(2))] 

\t<s<t+h / 

+ 1 - $ (A_(7)) + e^'^^^-") [1 - $ (A+(7))] 
= 1 - $ {x-{z)) + e2(A+^-)(5(*)--) [1 _ $ (x+(2))] + o{h), (22) 

where x±(^) := (git) - z)/^-^/^ ± (A + K-)h^/^. Returning to and using ([I2D, (HHD 
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and ( l22l) . we have 

r9{t) r _ n 

l2< / 2A,ig,t)igit)-z) 1 - $ + e2(A+™*)-) (1 _ $ (;^^(^))) dz 

Jg(t)-hi/4 L J 

+ 0{h) 

rait) 

=:2Ah{g,t) / {Ji{z) + J2{z))dz + o{h). 

Jgit)-h^/^ 

Consider the first term of the integrand (the second is treated in the same manner, since 
the exponential factor is (1 + o(l))). Denoting {K~ + fi) by K and changing variables 
y = {g{t) — z)/\/h, we have 



f9it) rh-^'' , _ . 

\ Ji{z)dz = h y [1 - <^{y - KVh)) dy 



'ait)~hUi 

<h I y [l-<^{y -KVh)jdy 



oo ^ 







/ {u + KVh){l-^u))du+ _ {u + KVh){l- <^{u))d 

Jo J-Ky/h 

< h Qe(Z2;Z > 0) +KVhE{Z;Z> 0) + K^h^ = ^h + o{h), 

where Z is a standard normal random variable. One can similarly show that the contri- 
bution from J2{z) will be of the same form. This gives 

h < Ah{g,t)h + o{h). 

Combining the bound with ffTTl) and (^lij, we conclude that 

h-'F,{T e{t,t + h)) < Ah{g, t) + o(l). 

Passing to the limit as /i ^ in the above expression yields a bound for the density of 
r, with the limiting value of Ah{g,t) given by B{g,t) defined in (fTUj) . This completes the 
proof of the theorem. □ 

For diffusions on R we can extend this result, considering the crossing of a lower boundary. 
Define 

g{t) := min 

— 0<s<t 

Then for the random variable p := inf {t > : Xt < g{t)} we have the following assertion. 
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Corollary 3.7. In case [A], let Xt and g(t) he as in Theorem \3.1l except with g{0) < x. 
Set 

L,:= inf {^^'(y)+^,^(y)). 

g(t.)<y<oo 

Then p has a density Pp{t) which satisfies 

Pp{t)<C{g,t), 

where 

C{g,t) := -{x - git) + x, (7(t))e^(^W)-^(^-)-(*/2)^W . 

t 



This assertion is a direct consequence of Theorem 13 .11 applied to the diffusion —Xt crossing 
the upper boundary —g(t). 

Similar ideas can be used to obtain lower bounds for the density of the first passage time, 
both for upper and lower boundaries, although again we are restricted to the case [A]. 
We use the same notation as previously, however instead of M{t) we define M*{t) by 

M*it):= sup {fi'iy) + fi'iy)). 

-oo<y<g{t) 

Theorem 3.8. In case [A], let Xt and g{t) be as in Theorem \3.1[ Assume that 

X < min {g{t) - K+t) . (23) 

Then the first passage time r of the process Xt of the boundary g{t) has a density Pr{t) 
which is bounded from below by 

p,it)>B%g,t), (24) 

where B*{g,t) is given by 

B*{g,t) := -{g{t) - K+t - x)g(t, x, ^7(t))e^(^«)-^(^)-(*/2)^*(*). 

Proof. The scheme of proof is the same as for Theorem 13.11 with straightforward changes 
to obtain a lower bound. In particular, instead of considering the linear boundary crossing 
for (7f(s), we consider similar probabilities for the linear function 

gl{s) = g{t) + K+{s-t). 

Inequality arises from the need for this boundary to initially be above x. □ 
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A similar result holds for lower boundaries as well. 



Corollary 3.9. In case [A], let Xt and g{t) be as in Theorem \3.8[ except with g{0) < x. 
Define L*{t) by 

V{t):= inf {^i'{y) + Ay))- 

g{t)<y<oo 

If X > supo<«i {g{t) + K^t) , then p has a density Pp(t), which is bounded from below by 

Pp{t)>C*{g,t){g,t), 

where 

C*{g, t){g, t) ■=lix- g{t) + K+t)q{t, x, g{t))eGi9(t))-G(x)-mwm) _ 



Similarly to the case of Corollary 13.71 the above assertion is a direct consequence of 
Theorem 13. 8[ considering the diffusion —Xt crossing the upper boundary —g{t). 

Remark 3.10. As we said in Section O the choice of the reference processes used is largely 
arbitrary within some constraints. If we have further information for a certain diffusion, 
the above results could be improved by substituting this diffusion as the reference process. 
In particular, several processes have the same bridges as either Brownian motion or Bessel 
processes, see |3j and references therein. In other words, if Xg is a process with the same 
bridge as a Bessel process, we have 



P„. ( sup (X, -gt{s))<0 

•.0<s<t 



X, 



P.. 



sup {Rs - gt{s)) < 

0<s<i 



Rf 



This could eliminate the need to use Lemma 13.31 proceeding straight to (fT6|) or its equiv- 
alent Lemma 13. 5[ This eliminates the exponential term in (|TU|) , however it also relies on 
knowing the transition density for the diffusion in question. 

Remark 3.11. Further, it is bounds for these conditional probabilities that limit the results 
obtained above and their accuracy. In particular. Lemma gives an excessive bound, of 
which tightening could improve on the bounds presented here. In case [B] it is also this 
lemma that needs to be extended to obtain a lower bound on the crossing density and to 
obtain results for lower boundaries. 



Observe that the above results can be extended to some time-inhomogeneous diffusions as 
well. In general, the transformation leading to a constant diffusion coefficient can result 
in a drift term which inhibits the use of the above methods. Even assuming we begin 
with a process with constant diffusion, a similar problem is encountered with the Girsanov 
density in Lemma [3.31 

However, for some time-inhomogeneous diffusions there is an obvious approach to avoid 
these complications. If our time-inhomogeneous process Yt is a function of a time- 
homogeneous diffusion Xt and time t, the idea is to manipulate the boundary crossing 
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probability so as to incorporate the time-inhomogeneity into the boundary. As a simple 
example, if Yt = r{t)Xt, with r{t) > 0, then we can write 

( sup (Ft - git)) > = f sup fx, - ^] > o] 

\0<t<l J \o<t<i V ^(^)/ / 

and proceed as before. 



4 Approximation Rates for Boundary Crossing 
Probabilities 



In this section we use the bounds of the previous section to extend the results of [5] from 
the Brownian motion case to general diffusions. The problem considered is calculating the 
probability that a given diffusion will stay between two boundaries over the time interval 
[0, 1]. That is, for a diffusion Xf and two functions g± on [0, 1] such that g-{0) < x < 
g+{Q)i we need to evaluate 

P, {g^{t) <Xt< g+{t),te [0, 1]) =: P{g^,g+) . (25) 

Again, this can be extended to the time interval [0,T], by a change of scale. We use the 
same notation for one-sided boundaries by replacing the redundant boundary with ±oo. 

An exact calculation of the above probability is possible only for a small number of 
diffusions and boundaries. Therefore in the general case, it is important to have an upper 
bound for the difference between probability (1231) and the corresponding probability for 
a pair of simpler boundaries which approximate the original ones and, at the same time, 
make computing the boundary crossing probability feasible. 

We denote the uniform norm of a function on [0, 1] by || ■ ||, that is \\g\\ = supo<t<i \git)\. 

Define De{g) to be the probability that the process Xt crosses the upper boundary g at 
some point t G [0, 1] without crossing the boundary + e, e > 0, on that time interval: 

D,{g) ■.= P{-oo,g + e)-P{-oc,g). 

Define D~ (g) similarly for lower boundaries: 

D;{g) := Pig - e,oo)- Pig, ^). 

Theorem 4.1. Let g± each satisfy condition ((H]), and consider functions f± on [0, 1] such 
that \ \g± — /ill < e for some e > 0. For to e [0, 1), let 

Btoig) ■= sup Big,t) and Ct^ig) := sup Cig,t), 

to<t<l to<t<l 
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yU := inf fi{x) and := sup fJ'i^), (26) 

— e<x<g{0)+K++e g{0)~K- -£<x<oo 

and set K* := max{0, — yu} and K := max{0, K" + Jt}. If fi > —oo, then 
\P (-00, g+) - P (-00, /+)| < max {D,ig+), D,{g+ - e)} , 

where 

De{9) < (^Tj^ + 2BMVT^^ + 2/r j e. (27) 

In case [A], if we also have JI < 00, then 

\P {g., g+) -PU-, /+) I < max {D,(^+) + D;ig_), D,{g+ - e) + D;{g. + e)} , 



where 



2 / 1 



D:{g) <\\l- y-^f^^ + 2CtMVT^, j+2Kje. (28) 

Remark 4.2. (i) Note that the two boundaries g± may satisfy condition ([8]) for different 
constants K^. For the purposes of this theorem, we have assumed the constants are the 
same, but if further accuracy is required, this assumption can easily be dropped. 

(ii) The choice of to is arbitrary and can be used to minimise the upper bounds for specific 
diffusion processes and boundaries. 

(iii) By setting 

B{g) := sup B{g,t) 

0<t<l 

we obtain the simpler, but less accurate, bound 



TT 



and similarly for {g). 

Proof. We begin by using an argument similar to the one employed in the proof of The- 
orem 1 in Since \ \g± — /±|| < we have 

P{g^+6,g+-e)<P{f^,f+)<Pig^~e,g+ + e). (29) 
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Further, 

P{g^,g+) < P{g--e,g+ + e), 

and hence 

< [P {g- -6,g+ + e)-Pig_,g+ + 6)] + [P (g^ , g+ + e) - Pig., g+)] . 
Each of these terms can now be bounded. For the second term, 



P{g.,g+ + e)-P{g.,g^)=F40< sup {Xt - g^t)) < e M (Xt - g^t)) > 

o<t<i o<t<i 



<PjO< sup {Xt-g+{t))<e 
V o<t<i 

= P{-oo,g+ + e) - P{-oo,g+) = D,{g+), 

while for the first term 

P{g~ -e,g+ + e) - P{g^,g+ + e) 

= fJ sup {Xt - g+{t) - e) < 0, < inf {Xt - g-{t)) < 
\o<t<i o<t<i 



<P. ^-e<mf^^{Xt-g-{t))<0 
= P{g-- 6, oo) - P{g^, oo) = D;{g^). 

Hence we have 

0<Pig^-e,g+ + 6)- P{g-,g+) < D,{g+) + D~{g.). (30) 

Using the same argument gives 

Q<P{g.,g+)-P{g.+s,g+-e)< D,{g+ ~ s) + {g. + e). (31) 

Now from (I29l)-(l3l]) we have 

\P{g.,g+) - PU'~, f+)\ < max {D,{g+) + D;ig^), D,{g+ - e) + D'ig^ + e)] . 

A simpler argument also gives 

|P(-oo,(?+)-P(-oo,/+)| < max {/},((?+), /},((?+-£)}. 

To complete the proof of the theorem, it remains to establish ( 1271) and ( l28l) . Since these 
bounds are derived in the same way, we will only prove the former one. 
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Note that 



De{g) = Px. < sup {Xt - git)) < e 
V o<t<i 

= f P.(r e dt)P, { sup (X, - g{s)) < e\Xt = g{t) 

Jo \t<s<l 

< [ P,(r G dt)F ( sup (X, - K+{s - t)) < e + g{t)\Xt = g{t)] 

Jo \t<s<l J 

= [ P,(r e cit)P,(t) ( sup (X, - K+s) < e + g{t)] . 

Jo \0<s<l-t J 



(32) 



Consider the integrand in the last hne. If we have a process Ys such that < Xg a.s. 
prior to the time inf{s > : > g{t) + e + K^s}, then 

Pg(t) ( sup (X, - K+s) <e + git)] < Pg(t) ( sup (F, - K+s) <e + git)] . 

\0<s<l~t J \0<s<\-t J 

To construct a process Yg meeting the above requirement, simply put 

dYs = fids + dWs, s > 0, Yq = git), 
with fi defined by fpUl) . For the new process, we clearly have 

P,(t)f sup (Ys-K^s) <e + git)] =f( snp {Ws - iK+ - fi)s) < e] 

\0<s<l-t J \0<s<\-t J 

< P ( sup iWs - K*s) < e] . 

\0<s<l-t / 

From here we continue as per the derivation of relation (14) in [5]. As is well-known (see 
e.g. (1.1.4) on p.250 in |4]), 

P( sup iWs-K*s)<e 

\0<s<l-t 



< $ ( K*VT^ + - ^ [K*VT^ - ^= ) + (1 - e-2^*^) 



Since <l>'(x) < i27i)-^^^ and 1 - e'^^*^ < 2K*e, we conclude that 

\0<s<l-t 
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( sup (X, - K+s) < e + git)] < J e + 2K*e, 

\0<s<l-t J V TTl^i — t) 



and therefore (1321) implies that 

Fix to £ [0, 1) and break the integral in ( l33l) as follows: 



TT V\/l -^0 Jto V^-t 



Using the bound ([9]) for the density Pa.(r G dt)/dt, with Btg{g) replacing B{g,t), then 
gives 

2 / 1 



+ 25i„((?)v/r^ +2K*. 



The derivation of (1281) uses the same argument. Theorem 14. II is proved. □ 
Corollary 4.3. For a Borel set B, define 

P{g.,g+-B):= {g^{t) < < t G [0, 1]; Xi G 5) . 

Then, under the assumptions of Theorem \A.l[ we have in the general case 

\P (-00, g+; B)-P (-cx), /+; B) \ < max {D,{g+), D,{g+ - e)] , 
while in case [A] 

\P ig^,g+- B)-P (/_,/+; B)\ < max {D,{g+) + D;{g^), D,{g+ - e) + D;{g_ +6)}, 
where D,. and satisfy the bounds given in Theorem 14. 1[ 

Proof. Essentially the same argument as for Theorem 14.11 gives the desired result. □ 
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As a simple consequence, we have an analogue to Corollary 1 of [5], which shows that 
using approximation by piecewise-linear functions results in a rate of convergence that is 
a quadratic function of the partition rank. 

Corollary 4.4. Let g± be continuously differentiable on [0,1], and g'_^ be absolutely con- 
tinuous, satisfying l^?^] < ^ < cxd almost everywhere. If = t^ < ti ■ ■ ■ < tn = 1 is a 
partition of [0, 1] of rank 6 = maxo<j<„ \ti — ti-i\, and f± are piecewise linear with nodes 
at the points (ti,g±(ti)), then 

\P(-oc,g+)^P{-oo,f+)\<^^5'. 

In case [A] we also have 

D 



\P{g.,g^)-P{f.J^)\<j^6\ 

where C and D can be bounded explicitly using Theorem \A.l[ In particular, for a uniform 
partition with ti = i/n, < i < n, the convergence rate of the above probabilities is 

Proof. Using Theorem I4.H we need only to show that 

e = \\f-g\\<l^5\ 

the proof of which can be found on p. 90 of [5j. □ 



Note that a similar result holds for the probabilities P{f-, f-\B). 



5 Examples 

In this section we illustrate the precision of the results from the previous two sections. 
Our bounds from Theorems 13.11 and 13.81 for the first passage time density are compared to 
some known densities for the Brownian motion, Bessel processes and Ornstein-Uhlenbeck 
processes. We also illustrate by numerical examples the results of Section H] and, and in 
the case of the Brownian motion, compare these with known results from the literature. 

5.1 Brownian Motion 

First we consider a Brownian motion starting at zero. Observe that for a linear boundary 
the bound of Theorem 13.11 is an exact result. Next we consider a curvilinear Daniels 
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boundary (as it has a known first passage time density) given by 

t , I Ki 



-45^ /t 



(34) 



where 5 7^ 0, ki > and /t2 G M subject to /tf + 4^2 > 0. For this boundary, the first 
passage time density is given by 



(see [7j for further information). We choose for our example the parameter values 6 = 
Ki = K2 = 0.5. Figure 15.11 displays graphs of the boundary (IMl) and the first passage 
time density ( |35ll compared to the upper and lower bounds as calculated by and ([2] 
respectively. 



The Daniels Boundary 



True density and its upper and iower bounds 
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Figure 5.1: The Daniels boundary with 5 = ki = K2 = 0.5, alongside its true first passage 

time density and upper and lower bounds. 

We can approximate the Daniels boundary by a linear function on [0, 1] with a uniform 
error of approximately 0.05, and clearly a piecewise linear boundary will reduce this. The 
results from Section H] for an approximating function / with || / — (?!! < s give the bound 

\P {-00, g)-P {-00 J)\<3.05e 

with to ^ 0.555. By way of comparison, the results in |5_j bound the above difference by 
3.73e. 

Note that choosing different parameter values to make the boundary flatter can increase 
the accuracy of the result. For example, with the choice 6 = 1.0, ki = K2 = 0.9 we have a 
maximum relative error for the density upper bound of less than 2.9% and the difference 
of the boundary crossing probabilities is bounded by 1.7e. 
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5.2 The Bessel Process of Dimension 4 



For a Bessel process of dimension 4, we choose the straight hne boundary g{t) = 1, and a 
starting point x = 0.5. The graph of the true first passage time density (given by (2.0.2) 
in [1]), as well as the upper bound calculated using Section [3l is shown on a semi-log scale 
in Figure \^72[ 




0.2 0.4 0.6 0.8 1 

Time 



Figure 5.2: First passage time density for a Bessel process and a straight line boundary, 
alongside its upper bound on a semi-log scale. 

As mentioned, improvements in Lemma [331 will help to improve this result, which anyway 
gives a useful bound on the error in the boundary crossing probabilities. Using a value of 
to ^ 0.75 gives the bound, for an approximating function /, with ||/ — (^H < £ < 0.05, 

\P{-oo,g)-P{-ooJ)\<2.25e. 

5.3 The Ornstein-Uhlenbeck Process 

Now consider an Ornstein-Uhlenbeck process St, which satisfies the SDE 

dSt = -e{St - I2)dt + dWt, 
where 6 > 0, fi E M.. According to [6J, for the 'hyperbolic' boundary 

g{t) = fi + Ae-^' + Be^\ 
where A,BeM. are such that x < fi + A + B, the first passage time density is given by 

Prit) = 2^^^t#^^p(t, X, g{t)), (36) 
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where p{t, x, z) is the transition density of the process St- 

^((2-/i)e^*-(x-/i))'\ 
1 - e^e* j ■ 

For the parameter choice 9 = 0.5, fi = 2,x = l, A = l and B = —1, Figure ESI displays 
the true first passage time density ( !36l) compared to the upper bound calculated by ([9]). 



p{t,x,z) 



7r(e2^* - 1) 



1/2 



exp 



4 




C 

CD 
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Figure 5.3: First passage time density of a hyperbolic boundary for an Ornstein-Uhlenbeck 

process, alongside its upper bound. 

Although this process is an example of case [B], we are unable to obtain a meaningful 
lower bound as M*{t) is unbounded. 

Using the results of Section H] with to ~ 0.28 gives the bound 

|P(-oo,(?)-P(-oo,/)|< 1.885, 

when ||/-^|| <e < 0.05. 
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